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Abstract 

The Gonihedric Ising model is a particular case of the class of models defined by 
Savvidy and Wegner intended as discrete versions of string theories on cubic lattices. 
In this paper we perform a high statistics analysis of the phase transition exhibited 
by the 3d Gonihedric Ising model with k = in the light of a set of recently 
stated scaling laws applicable to first order phase transitions with fixed boundary 
conditions. Even though qualitative evidence was presented in a previous paper to 
support the existence of a first order phase transition at k = 0, only now are we 
capable of pinpointing the transition inverse temperature at /3 C = 0.54757(63) and 
of checking the scaling of standard observables. 
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1 Introduction 



In a recent paper [1] we have studied the effects of freezing the boundaries in a 
Monte Carlo simulation near a first order phase transition. More specifically, 
we checked (and postulated one of) the scaling laws governing the critical 
regime of the transition by means of a Monte Carlo simulation of the 2d, 8- 
state spin Potts model. These new scaling laws, theoretically analyzed by C. 
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Borgs and R. Kotecky and by I. Medved [2], imply a major change in the 
critical behavior analysis. 

The MC simulation of a system with fixed boundary conditions (F.B.C.) in- 
stead of the standard periodic ones (P.B.C.) is more than a simple academic 
exercise. Indeed, the numerical analysis of the 3d Gonihedric Ising model re- 
quires fixing the spins of some internal planes. If periodic boundary conditions 
are adopted, the fixing of these internal planes is just equivalent to the simu- 
lation of the system in a box with fixed boundary conditions. For this reason, 
the Gonihedric Ising model with k — 0, which manifests a first order phase 
transition [3], needs to be reanalyzed in the light of the appropriate scaling 
laws. Moreover, in our recent paper [1], the new scaling laws were checked 
for a two dimensional system, so the 3d Gonihedric Ising model offers the 
opportunity to extend their verification to 3d lattices. 

In the present paper we perform a high statistics study of the 3d Gonihedric 
Ising model with k — at the transition point on lattices up to 20 3 . Our 
analysis of the scaling behavior of some standard thermodynamical magni- 
tudes (specific heat, susceptibility and energetic Binder cumulant) confirms 
the above-mentioned scaling laws and shows the importance of applying the 
correct scaling forms when fixed boundary conditions are present. 

This letter is divided as follows. A brief summary of the Gonihedric Ising 
model is contained in Sec. 2. The scaling laws for first order phase transitions 
are stated in Sec. 3, comparing the laws for fixed boundary conditions with 
their periodic counterparts. Sec. 4 and Sec. 5 are devoted to the numerical 
simulation and analysis of results and Sec. 6 summarizes the conclusions of 
our work. 



2 The Gonihedric Ising model at k — 

Adding extended range interactions, particularly with different sign couplings, 
to the standard Ising model in two and three dimensions gives a very rich [4] 
phase structure. One particular class of models with such extended interac- 
tions, the so-called Gonihedric Ising models, have recently aroused interest 
because of their putative connection with random surface models and strings. 
The original discretized random surface model was developed by Savvidy et 
al. [5] with the action 




<ij> 



(1) 
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where the sum is over the edges of some triangulated surface, 0((Xij) = |vr — 
CKjjl^, C is some exponent, and is the dihedral angle between neighbouring 
triangles with common link < ij >. It was christened the Gonihedric string 
model. 

The above action was translated to plaquette surfaces by Savvidy and Weg- 
ner [6,7] who rewrote the resulting theory as a generalized Ising model by us- 
ing the geometrical spin cluster boundaries to define the plaquette surfaces. In 
view of its relation to the Gonihedric string model, this new action was named 
the Gonihedric Ising model. In what follows we shall consider the three dimen- 
sional version of this model, whose Hamiltonian contains nearest neighbour 
(< i,j >), next to nearest neighbour (<< i,j >>) and round a plaquette 
([i,j,k,l]) terms 

H = 2kJ2 ~~ 2 + H WjWi- (2) 

<ij> Z Z [i,j,k,l] 

For generic couplings the spin clusters in the above Hamiltonian generate a 
gas of surfaces with energy contributions from area, extrinsic curvature and 
self- intersections [8]. A noteworthy feature of the particular ratio of couplings 
in Eq. (2) is the flip symmetry which is not present in the generic case. It is 
possible to flip any plane of spins at zero energy cost when T = 0, so the zero 
temperature ground state is degenerate, with any layered configuration being 
equivalent to the ferromagnetic state. A low temperature expansion shows that 
this symmetry is lost when T^O and k ^ [7]. k — however constitutes a 
special case - the flip symmetry remains even at finite temperature. 

There is agreement on the phase structure of the Hamiltonian in Eq. (2) 
from both Monte Carlo simulations and cluster- variational (CVPAM) meth- 
ods: when k > there is a single continuous transition from a paramag- 
netic high temperature phase to (with appropriate boundary conditions in 
the Monte Carlo case) a ferromagnetic phase. The simulations of Ref. [9] used 
fixed boundary conditions in order to define a magnetic order parameter; the 
reason was that it was found that with the use of standard periodic bound- 
ary conditions flipped spin layers, with arbitrary interlayer spacings, made it 
unfeasible. 

The nature of the transition for k ~ was then investigated in Ref. [3]. A 
zero temperature analysis [9] shows that there is a further "antiferromagnetic" 
symmetry in the ground state when k — 0, which is already apparent from the 
Hamiltonian itself. This extra symmetry, and the persistence of flip symmetries 
at non-zero T suggest that k = is a special point in the space of Hamiltonians 
Eq.(2). Even though the results of Ref. [3] suggested the presence of a first 
order phase transition at k — 0, a complete finite size analysis of the transition 
was not performed at that time for want of a better knowledge of the scaling 
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laws applicable with fixed boundary conditions. 



3 The new scaling laws for frozen boundaries 

As mentioned in the introduction, the scaling laws applicable to systems sim- 
ulated with fixed boundary conditions were deduced and studied in Ref. [1,2]. 
The numerical analysis of Ref. [1] was performed on the 2d 8-state Potts 
model. Since the difference between the corresponding scaling laws for fixed 
and periodic boundary conditions are highly volume-dependent, in addition 
to its intrinsic interest the simulation of the 3d Gonihedric Ising model is a 
good testing ground for the new scaling laws on a 3d lattice. 

A main feature of the F.B.C. simulations is the shift of the infinite volume 
inverse temperature by a 1/L correction term, caused by surface effects, in- 
stead of the 1/L d correction term due to volume effects seen in the periodic 
case. The same change in the shift is also observed for the energetic Binder 
parameter with fixed boundary conditions. 

Moreover, the surface corrections to the volume scaling of the specific heat 
and the susceptibility become of order L d ~ x in the fixed case instead of the 
almost negligible 1/L d . 

Table 1 summarizes the scaling laws for a first order phase transition for both 
periodic and fixed boundary conditions. 



4 Numerical simulation 

As we have already noted, the flip symmetry poses something of a problem 
when carrying out simulations since it means that a simple ferromagnetic order 
parameter 



will be zero, because of the observed layered nature of the ordered state. 
Staggered magnetizations are of no use since the inter layer spacing can be 
arbitrary. On a finite lattice it is possible, however, to force the model into 
the ferromagnetic ground state by fixing sufficient perpendicular spin planes, 
either internally if P.B.C. are used or on the boundaries of the lattice: both 
possibilities being exactly equivalent. 




(3) 
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As in our previous work [3], we choose to fix internal planes of spins in the 
lattice, while retaining the periodic boundary conditions. This has the desired 
effect of picking out the ferromagnetic ground state. We can therefore still 
employ the simple order parameter of Eq. (3). For k = the Hamiltonian we 
simulate is 1 

H = ^ WjOkOi- (4) 

[i,j,k,l] 



Table 2 summarizes the details of the simulations that have been performed 
from L = 10 up to L = 20. The lattice updating used a simple Metropo- 
lis algorithm. The number of production Monte Carlo sweeps varies from 
n prod = 20 000 000 for L = 10, to n prod = 200 000 000 for L = 20. We 
took measurements of the energy and the magnetization only every nm P = 4 
or ft-flip — 8 sweeps, and, consequently, the number of total measurements per 
run is n meas = n pro d/nfji P . We left at least 21nfli P r e thermalization sweeps 
before taking measurements [10]. To estimate the autocorrelation time of 
energy measurements r e , we use the fact that r e enters the error estimate 
ejK = ^2 r e /n meas e na i ve for the mean energy < E > of n mcas correlated en- 
ergy measurements of variance 

4ive = E (< E > -Ej) 2 /(n meas - 1). (5) 
j'=i 



The "true" error estimate ejK is obtained splitting the energy time-series into 
50 bins, which were in their turn jackknived [11] to decrease the bias in the 
analysis. 

In Fig. 1 we present the energy time-series for the L = 20 and @mc — 0.5064 
simulation run. The expected characteristic behaviour of a first order phase 
transition can be clearly seen. The system remains in one of the two coexisting 
phases for a long period of time. The energy histogram for the full series is also 
presented in the Figure. The similar height of the two peaks confirms that the 
simulation was performed very near the pseudo-critical inverse temperature. 

In addition to the qualitative analysis of the histograms, we have computed 
the specific heat, magnetic susceptibility and the energetic Binder parameter 
at nearby values of (3mc by means of standard reweighting techniques [12]. 
These observables are defined as 



It is perhaps worth emphasizing that spins live on the vertices of the cubic lattice 
rather than on the links, so the model of Eq. (4) is not the three dimensional Z2 
gauge model that is dual to the three dimensional Ising model. 
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C{(3) = ^{(E*)-(E)\ (6) 
X(/5) = ^((M 2 )-(M) 2 ), (7) 

m^-^yr (8) 

In Table 3 we show the extrema of the magnitudes defined above, together with 
their pseudo-critical inverse temperatures. The error bars of these quantities 
have been estimated splitting the time-series data into 50 bins, which were 
then jackknived to decrease the bias in the analysis of reweighted data. 



5 Analysis of results 



Once we have the results from the numerical simulation on finite lattices, we 
can proceed to analyze the data by fitting to the scaling laws of Table 1. 

In Table 6 we show the results of fitting the pseudo-critical /3s of C max , Xmax 
and B min to the ansatz 

/WL)=/3 c + | + |f (9) 



suggested by the finite-size scaling laws presented in Table 1. For x max and 
_B mm the fits were rather poor if L = 10 was included, so it was discarded. 
For C max both sets L = 10 — 20 and L = 12 — 20 were fitted. Focusing on the 
L = 12 — 20 fits, we can discern only very minor differences in the estimated 
P c depending on the observable used to extract it. These are so small that we 
can safely average to obtain 

(3 C = 0.54757 ±0.00063 (10) 



Since the /3 c 's extracted from the three observables were not independent, we 
have kept the error bar common to them all. In Fig. 2 we depict the fit for 
^(1) in the range L = 10 — 20. The error bars in the Figure are so small 
that they show up only as horizontal dashes. 

The results of the fits to the specific heat and susceptibility maxima, C max and 
Xmax, together with the energetic Binder parameter minimum are summarized 
in Table 6. The goodness-of-fit, Q, is excellent for the three observables. 

Note that the surface correction coefficients ai and b\ are, in absolute value, 
from one to two orders of magnitude larger than the coefficients a 2 and b 2 of 
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the dominant contribution V = L 3 . It is precisely this fact which makes it 
necessary to use the scaling ansatz C max (L) = a + a\ L? + a<i L 3 , and allows 
us to estimate the corrections to the leading term. 

6 Conclusions 

We have performed a numerical simulation of the 3d Gonihedric Ising model 
at k = in order to determine the thermodynamic characteristics of its phase 
transition. Previous analysis suggested the existence of a first order phase 
transition, but a complete finite size analysis of the transition was not carried 
out. The special features of this model, which requires a simulation where 
three perpendicular spin planes need to be fixed during the simulation, do not 
allow a direct application of the standard finite size scaling laws for periodic 
boundary conditions at a first order transition. In fact, to keep these planes 
fixed is equivalent to performing a simulation with fixed boundary conditions 
(F.B.C.), giving rise to the need for a different set of scaling laws. They were 
reviewed in Sec. 3. Our numerical analysis of the thermodynamic quantities 
has shown that the critical behavior of the 3d Gonihedric Ising model is per- 
fectly described in terms of F.B.C. scaling laws. As a result of this work, we 
have been able to accurately determine the inverse critical temperature of the 
model, i.e. f3 c = 0.54757(63). Furthermore, our simulation has extended the 
verification of the F.B.C. scaling laws to a three dimensional lattice model. 
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Table 1 

Scaling laws for Periodic versus Fixed Boundary Conditions. 





P.B.C. 


F.B.C. 


jf c eaks (L) = 

Cmax {L) = 
Xmax(L) = 
Brnin(L^) = 


/3 c (oo) + ^ + 0(^) 
70 + ^L d + 0(±) 
5 + S 2 L d + 0( 

^o + fi + 0(^) 


/3 c (oo) + ^ + 0(^) 
co + csL^ + O^- 1 ) 
eo + esL^ + O^- 1 ) 



Table 2 

Monte Carlo parameters of the simulation. L 3 is the lattice size, n t i lCTm the number 
of Monte Carlo sweeps during thermalization (in thousands), and n pro d the number 
of production runs (in millions). Measurements were taken every n^p = 4 Monte 
Carlo sweeps for all the simulations, except the latest; the number of bins was 50. 



L 



PMC "-therm "-prod "flip 



"thcrm/"-flip "prod/"flip 
To 2 T e 



10 



12 



14 



15 



18 



20 



0.4580 



0.4748 



0.4864 



500 



500 



500 



20 



20 



20 



25 



45 



278 



0.4910 500 20 4 1 011 



0.5013 2 500 22 4 24 871 



0.5064 36 500 200 8 216 098 



5 000 100 000 



2 778 



450 



25 



21 



55 556 



8 993 



124 2 473 



111 



58 



Table 3 

Extrema for the (finite lattice) specific heat, C max , the susceptibility, Xmax? and 
the energetic Binder parameter, B m \ n , together with their respective pseudo-critical 
inverse temperatures. 



L 


C 


Cmax 




Xmax 


A^mm 


-^min 


10 


0.457919(21) 


5.6945(79) 


0.456842(22) 


7.042(12) 


0.455064(21) 


0.638537(47) 


12 


0.474753(16) 


12.120(21) 


0.474470(15) 


16.305(31) 


0.473468(16) 


0.635656(61) 


14 


0.486349(21) 


25.172(45) 


0.486275(21) 


36.264(73) 


0.485647(21) 


0.628430(85) 


15 


0.490922(30) 


34.900(76) 


0.490884(30) 


51.91(13) 


0.490374(30) 


0.62432(12) 


18 


0.501280(72) 


78.99(38) 


0.501273(72) 


128.10(69) 


0.500979(72) 


0.61246(36) 


20 


0.506366(69) 


121.57(52) 


0.506364(69) 


206.61(99) 


0.506149(69) 


0.60620(36) 



10 



Table 4. Pseudo-critical inverse temperature fits. Q is the goodness- of- fit. 



range L's 


P^ x (L)=f3 c + ai /L + a 2 /L 2 


= (3 c + d 1 /L + d 2 /L 2 


^ ax (L)=/3 c + ei /L + e 2 /L 2 




Q 0c ai a 2 


Q Pc ai a-2 


Q Pc ai a 2 


10 - 20 


0.89 0.54868(34) -0.7848(84) -1.229(51) 






12 - 20 


0.73 0.54867(63) -0.785(18) -1.23(12) 


0.85 0.54730(63) -0.736(18) -1.65(12) 


0.86 0.54674(63) -0.711(18) -2.02(12) 



Table 5. Fits on the extrema of C max , Xmax and B m { 



range L's 


C m ax(i) = a + a 1 L 2 + a 2 L 3 


Xmax(^) =b + b 1 L 2 + b 2 L 3 


B ulin (L) = B + B 1 /L + B 2 /L 2 






Q b h b 2 


Q Bq Bi B 2 


10 - 20 
12 - 20 


0.16 14.43(17) -0.4434(36) 0.03561(20) 
0.098 14.79(52) -0.4491(83) 0.03587(40) 


0.62 39.92(92) -1.035(15) 0.07257(73) 


0.014 0.4992(15) 2.851(35) -14.57(21) 
0.21 0.5065(30) 2.643(84) -13.12(57) 




1 2 

histogram(E) 

10 measurements 



Fig. 1. Energy time series and corresponding energy histogram for L = 20 and 
Pmc = 0.5064 
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0.4 



15 



25 



Fig. 2. Finite-size scaling analysis of the pseudo-critical /?^ ax in the range L = 10— 20 
by means of the ansatz /3^ ax (L) = (3 C + a±/L + (12/L 2 . The infinite volume critical 
point obtained from the fit is [3 C = 0.54868(34), with a goodness-of-fit Q = 0.89. 
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